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ABSTRACT 

We investigate the ability of spectroscopic techniques to yield realistic star formation histories 
(SFHs) for the bulges of spiral galaxies based on a comparison with their observed broadband colors. 
Full spectrum fitting to optical spectra indicates that recent (within ~ 1 Gyr) star formation activity 
can contribute significantly to the y-band flux, whilst accounting for only a minor fraction of the 
stellar mass budget which is made up primarily of old stars. Furthermore, recent implementations of 
stellar population (SP) models reveal that the inclusion of a more complete treatment of the thermally 
pulsating asymptotic giant branch (TP-AGB) phase to SP models greatly increases the NIR flux for 
SPs of ages 0.2-2 Gyr. Comparing the optical-NIR colors predicted from population synthesis fitting, 
using models which do not include all stages of the TP-AGB phase, to the observed colors reveals that 
observed optical-NIR colors are too red compared to the model predictions. However, when a 1 Gyr 
SP from models including a full treatment the TP-AGB phase is used, the observed and predicted 
colors are in good agreement. This has strong implications for the interpretation of stellar populations, 
dust content, and SFHs derived from colors alone. 

Subject headings: galaxies: bulges — galaxies: evolution — galaxies: formation — galaxies: stellar 
content 



L INTRODUCTION 

A detailed understanding of the stellar populations 
(SPs) that make up the integrated spectral energy dis- 
tributions (SEDs) of both local and distant galaxies can 
provide important constraints for models of galaxy for- 
mation. However, confronting observations with SP syn- 
thesis models for the purpose of translating the former 
into physical parameters, is not without significant chal- 
lenges. While SP modeling has seen tremendous progress 
over the past decade (e.g., Bruzual & Chariot 2003, here- 
after BC03; Le Borgne et al. 2004; Maraston 2005, here- 
after Mar05; Schiavon 2007), misinterpretations in the 
data-model comparison are still common. 

Broadband colors are often used as a proxy for SP pa- 
rameters. The well-known s.ge/Z degeneracy in optical 
colors is greatly lifted by the addition of NIR bands. 
However, in broadband-based analyses, extinction and 
reddening effects from interstellar dust cannot be ruled 
out. Indeed, observed optical-NIR colors that lie red- 
ward of the model grids, i.e., in a region not supported 
by a naked SP of any age/Z/SFH, have typically been 
attributed to dust reddening (e.g., Peletier et al. 1999; 
MacArthur et al. 2004; Carollo et al. 2007). Meanwhile, 
red optical-NIR colors within the model grids are as- 
sumed to result from old and metal-rich SPs. 

In light of recent implementations of SP models which 
account for a more detailed treatment of the effects of 
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the thermally pulsating asymptotic giant branch (TP- 
AGB) stellar evolutionary phase (e.g., Mar05; Coelho 
et al. 2007; Lee et al. 2009), a reassessment of the ob- 
served red optical-NIR colors in galaxies is in order. 
Accounting for the TP-AGB is accomplished in Mar05 
through empirical calibration to LMC globular clusters 
which have independent age and Z measurements. Alter- 
natively, synthetic model tracks can be used. The recent 
tracks of Marigo & Girardi (2007), whose predictions are 
in good agreement with those of Mar05, account for nine 
evolutionary stages (including the C-, M-type, and super- 
wind mass loss phases) in the TP-AGB. (Note that the 
older tracks used in the BC03 models account for only 
a single evolutionary stage for each evolutionary phase 
in the TP-AGB.) In these models, the TP-AGB phase 
is active in SPs of ages ~ 0.3-2 Gyr and contributes sig- 
nificantly to the NIR fiux (accounting for up to ~ 80% 
of the NIR flux), leading to very red optical-NIR colors 
for this age range. These models have already been used 
to resolve the uncomfortably large stellar masses and old 
ages derived for high-z galaxies using models which do 
not account for all phases in their treatment of the TP- 
AGB phase (Maraston et al. 2006). 

In the meantime, full spectrum fitting techniques have 
proven to be effective at recovering the underlying stel- 
lar content of integrated galaxy spectra (e.g.. Heavens 
et al. 2000; Cid Fernandes et al. 2005; Walcher et al. 
2006; MacArthur, Gonzalez, & Courteau 2009, hereafter 
Mac09). This type of analysis provides a stochastically- 
sampled SFH for a given integrated spectrum and allows 
for true average, as opposed to single SP-based, age and 
Z estimates. Full spectral fitting is typically limited to 
optical wavelengths, as set by the observations. However, 
as most models provide SEDs in the UV-FIR range, one 
can use such model fits to make predictions for the galaxy 
light in other wavelength regions. For example, a com- 
parison of predicted versus observed optical-NIR colors 
could provide valuable insight into both the reliability 
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of the spectral fits as well as any shortcomings of the 
SP models themselves. This is precisely the approach 
we carry out here using our analysis of Gemini/GMOS 
long-slit spectra of local spirals (Mac09) and new optical 
imaging for the same galaxies from the Palomar Obser- 
vatory and NIR imaging from the Two Micron All Sky 
Survey (hereafter 2MASS; Skrutskie et al. 2006). As we 
are interested in predictions from models that use dif- 
ferent treatments of the TP-AGB phase, we focus on 
two galaxies from our spectroscopic sample (NGC 628 & 
UGC 2124) with the highest quality data and whose de- 
rived SFHs indicate significant contribution (in y-band 
light-weight) from a 1 Gyr SP, i.e., where the TP-AGB 
NIR signature is expected to be strong. 

2. DATA 

The long-slit spectroscopic data used for this study, 
from Mac09, were collected with the Gemini Multi- 
Object Spectrograph (GMOS; Hook et al. 2004) on 
the 8-m Gemini North telescope. The GMOS detector 
and B600 grating combination provides a spatial resolu- 
tion of 0.072"/pix and a dispersion of 0.45A/pix. The 
slit field of view (FOV) of 5' x 2" provided a boxcar 
width = 10.81 A resolution. The spectral coverage spans 
-4050 -6750 A. 

Optical imaging in the BVRI Bessel filters was col- 
lected at the Palomar 200-inch telescope in 2006 De- 
cember. The Large Format Camera consists of a six 
CCD mosaic with a FOV of 6x12' per CCD. Our galax- 
ies have optical diameters < 10.5', leaving enough room 
beyond the optical radius on a single chip for sky mea- 
surements. Typical exposure times were 1-2 minutes and 
2x2 binning gave 0.363"/pixel. The seeing ranged from 
FWHM = 1.8-3.7". The data were bias-subtracted and 
flat-fielded using a combination of night-sky and dome 
flats to carefully account for both the high-frequency spa- 
tial sensitivity (dome) and the large-scale illumination 
(sky) of the chip. Due to the bright-moon conditions, 
it was necessary to model the background with a sur- 
face. This induced a typical sky error of ~ 0.01 mag. The 
data were calibrated following the technique of Courteau 
(1996) using stars from several Landolt fields (Landolt 
1992), for typical calibration errors of ^ 0.03 mag. 

Imaging in the NIR JHKs filters were obtained from 
the 2MASS online database. The background levels were 
carefully re-measured for each image and the photometric 
zeropoints were taken from the 2MASS headers. 

All data (imaging and spectroscopic) were corrected 
for Galactic foreground extinction (Schlegel, Finkbeiner, 
& Davis 1998). 

2.1. Radial Extractions 

In order to make a direct comparison between informa- 
tion gleaned from optical spectra versus broadband imag- 
ing, we must compare the spectral fits with colors from 
the same effective radial bins. As such, rather than the 
usual azimuthal extraction of SB profiles (e.g.,, Courteau 
1996; MacArthur, Courteau, & Holtzman 2003), profiles 
were extracted from all bands [BVRIJHKs) using the 
same position angle and slit width as for the GMOS ob- 
servations. We accounted for differences in the observed 
point spread functions by convolving all profiles (imaging 
and spectroscopic) with the largest FWHM of all obser- 



vations for each galaxy. Finally, radial binning for both 
the imaging and spectroscopy was set to the largest mea- 
sured dispersion, which is that of 2MASS (l"/pixel). 

3. AGE, METALLICITY, AND SFH FROM 
POPULATION SYNTHESIS 

The SP fitting technique used here is as described in 
Mac09; a brief summary is provided below. Our "full 
population synthesis" technique consists of a bound con- 
strained optimized fit representing the relative contribu- 
tion of each of the 70 model templates to the observed 
spectrum. The only fixed bound is that of no nega- 
tive template contributions. The templates are from the 
BC03 models which provide SEDs representing simple 
stellar populations (SSPs), i.e., single bursts of star for- 
mation (SF) at a given age and Z, at a resolution of 
- 3 A FWHlVfl. The selected library of 70 SSP templates 
covers the age range 0.001-20 Gyr and metallicities of 
Z = 0.0004-0.05 and we use the models with the Chabrier 
(2003) initial mass function. Our procedure also allows 
for dust reddening to the observed SED according to the 
prescription of Chariot & FaU (2000). See Mac09 for 
details. To accommodate any non-stellar contributions 
to the observed spectrum not considered in the BC03 
models, we used an iterative masking scheme whereby 
deviant | data-model] points are given zero weight in the 
fit. 

Figure [T] shows the full population synthesis fits to the 
central spectra of NGC 628 and UGC 2124 (extracted 
as described in Section [O]) . and in Figure [2] we plot the 
inferred SFHs from the fits. In Table [1] we list for each 
galaxy the percent-light and mass contributions of "very 
young" (0.001-0.4 Gyr), "young" (1-2 Gyr), "intermedi- 
ate" (4-7 Gyr), and "old" (10-20 Gyr) age SSPs to the 
fits. 

For both galaxies there is a contribution from 13 to 
20 Gyr SSPfl While their weight in light may not seem 
dominant, these old populations contribute 74% and 93% 
to the stellar mass budget of NGC 628 and UGC 2124, 
respectively (Table [T]). However, we are particularly in- 
terested here in a significant contribution from SSPs that 
would be expected by the Mar05 models to have NIR ex- 
cesses from their treatment of the TP-AGB stellar evo- 
lution phase. The age range for which the signature is 
present is highlighted by the gray shading in Figure [2] 
and referred to as "young" in Table [TJ The contribution 
to the y-band normalized light of the young population 
of interest is 67% and 33% for NGC 628 and UGC 2124, 
respectively. Thus, if our SFHs are reliable, according to 
the Mar05 models, these galaxies would be expected to 
show signs of the TP-AGB NIR excesfl 

4. "PREDICTED" VERSUS OBSERVED 
OPTICAL-NIR COLORS 

^ Although see §3.4.1 of Mac09 for a detailed analysis of the true 
resolution of these models. 

^ We are not concerned here with model ages that are older than 
the age of the universe as absolute model SSP ages are not precisely 
calibrated. 

* Note that we did not perform the population synthesis fits 
with the Mar05 models as they have a lower resolution than BC03, 
thus reducing the predictive power of the spectra to disentangle the 
high frequency features in the fitting process. A direct compari- 
son between fits using different models awaits the release of higher 
resolution models. 
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Fig. 1. — Central observed spectra (black) and full population synthesis fit (red) for NGC 628 (left) and UGC 2124 (right). Gray shading 
indicates regions masked in the fit as determined by our iterative "cr-clipping" procedure as well as the CCD gap regions (green vertical 
dash-dotted lines) which are always masked. Shown at lower right on each panel are the average light-weighted age, (A)j , metallicity, (Z); 
, effective Ty, and of the fit. The bottom panels show the percent data— model residuals. In the bottom panels, dashed vertical lines 
indicate variable sky-lines, and dotted vertical lines indicate emission lines prevalent in star forming (H II) regions. 



TABLE 1 

Percent Contributions of all SSPs in Given Age Ranges (Very Young = 0.001-0.4 Gyr; Young = 1-2 Gyr; 
Intermediate = 4-7 Gyr; Old = 10-20 Gyr) to Full Population Synthesis Fits Weighted by Light (V-band normalized) and 

Mass. 
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Fig. 2. — Light-weighted SFHs of the central spectra for NGC 
628 and UGC 2124 derived from the full population synthesis fits. 
The dotted horizontal lines guide the eye for each galaxy (labeled 
at right). The horizontal axis is the SSP age. Point size is propor- 
tional to the relative weight in the fit (normalized to the V-band) 
and the colors and point types code the SSP metallicity, as indi- 
cated. The gray shading indicates the age range within which the 
TP-AGB phase is active. 



The full population synthesis fits to the spectra are 
constrained by the observed 4050-6750 A region (gray 
shading in Figure [3]) and, modulo emission lines, repre- 
sent this region very well (see Figure [Ij. Since the BC03 
models cover the range 90 A -160 /zm, we can use the 
model fits to "predict" colors in other bands. Any dif- 
ference between colors predicted in this fashion and the 
observed colors can be explained by any of the following 
scenarios: unrealistic SFHs due to degeneracies in the 
age/Z/dust plane that cannot be constrained with op- 



tical data alone; inadequacy of the template library to 
represent the full coverage of SP age and Z in any given 
integrated spectrum; unrealistic estimate of the dust ex- 
tinction from the population synthesis fits; errors in the 
SP modeling predictions outside the optical range due to, 
for example, errors in the stellar evolutionary physics. 

To test for any discrepancies, we now compare the pre- 
dicted optical-NIR colors from the full population syn- 
thesis BC03 model fits to the observed colors. A di- 
rect estimate of the predicted NIR excess is complicated 
by a number of issues related to differences between the 
BC03 and Mar05 models. First, for a given age and 
Z , the model predictions are not identical, even within 
the optical limits of our spectra. Second, the SEDs pro- 
vided in each model do not cover precisely the same 
age/Z grid. To illustrate differences between the BC03 
and Mar05 model SEDs, we compare in Figure [3] the 
lGyr/Z = 0.05 SSP from BC03, which contributes sig- 
nificantly to both SFHs in Figure [51 to the SED closest 
in age and Z from the Mar05 models (1 Gyr/Z=0.04). 
While most of the difference is attributed to the differ- 
ent treatments of the TP-AGB phase between BC03 and 
Mar05, some of it could also be due to the overall larger 
contribution of AGB stars in the Mar05 models. For this 
comparison, the BC03 models were smoothed with a box- 
car of FWHM 20 A and resampled to A A = 20 A/pix to 
roughly match the resolution and sampling of the Mar05 
SEDs. 

Given the above caveats, we now look at the observed 
versus predicted optical-NIR colors for NGC 628 in Fig- 
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Fig. 3.— Comparison of SSP SEDs for the Mar05 and BC03 
models. The dashed red Une is the lGyr/Z = 0.05 SSP from 
the BC03 models. The solid blue line is the 1 Gyr/Z = 0.04 SSP 
from the Mar05 models. The relevant filter response curves for 
our observed colors are shown as dark gray curves in the up- 
per panel. The lower panel plots the model differences, where 
%diff=100*{BC03-Mar05)/[(BC03-|-Mar05)/2]. The gray shad- 
ing represents the optical range of our CMOS spectra within which 
the population fits are constrained. 



ure 31 In all color combinations, the model fits (red 
open triangles) predict colors that are too blue compared 
to the observed (black solid triangles) optical-NIR col- 
ors. The blue arrows represent the difference in color 
between the BC03 1 Gyr/Z = 0.05 SSP and the Mar05 
lGyr/Z = 0.04 SSP. Indeed, the blue vectors represent 
the data-model difference very well. For UGC 2124 in 
Figure [3 it appears that a steeper vector would be re- 
quired to precisely match the data-model difference. In 
Mac09 we demonstrated by a comparison of absorption- 
line indices that there is evidence of an enhanced [a/Fe] 
SP in the bulge region of UGC 2124 (see Figure 4 in 
Mac09), which could also affect the observed colors (only 
solar-scaled abundance ratio models are considered here). 
To assess whether abundance ratios could compensate 
for the extra steepness in the data-model difference, we 
also plot in Figure [5] a green vector which represents 
the effect of a super-solar abundance ratio star at a 
given Te///log(g)/Z from the synthetic stellar models 
of Coelho et al. (2005) (see the figure caption for de- 
tails). The Coelho models do not extend far enough into 
the NIR for an accurate X-band estimate. A combined 
effect of both TP-AGB treatment and super-solar [a/Fe] 
would be represented by the addition of the two vectors 
and resulting in a steepening that more closely accounts 
for the data/model offsets. 

As a further check on the model predictions, we plot 
in Figure m NIR-NIR colors for both galaxies. Again we 
see a difference between the observations and model pre- 
dictions, but not in the same sense as the optical-NIR 
colors. While the J — H color is predicted to be redder 
in the Mar05 models, the J — K and H — K colors are 
predicted to be bluer. As for the optical-NIR colors, the 
Mar05 models could account for these differences. The 
only mismatch is in the H — K versus J — H colors for 
UGC 2124, but this could very well be another manifes- 
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Fig. 4. — Color-color plots for NGC 628. The red open trian- 
gles are colors measured ("predicted") from the BC03 model full 
population synthesis fits to the observed central spectrum. The 
black solid triangles are the corresponding observed colors. Obser- 
vational errors (calibration and sky subtraction) are indicated at 
lower right. Overplotted are grids of BC03 models with exponen- 
tial SFHs (see MacArthur et al. 2004 for details). Solid and dashed 
lines are iso-Z and iso-age, respectively, with darker shading in- 
dicating higher Zs and older ages; see middle bottom panels for 
labels of average age and Z. Blue arrows: model prediction dif- 
ferences between the lGyr/Z = 0.05 BC03 SSP (red dashed curve 
in Figure [Sl and the Mar05 lGyr/Z = 0.04 SSP (solid blue curve 
in Figure[3|. The arrows point toward the Mar05 predictions (from 
their respective location for the BC03 models). 
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Fig. 5.— Same as Figure g] but for UGC 2124 and the 
green arrows show the effect of a non-solar abundance ra- 
tio SSP on the colors. Specifically, this represents the differ- 
ence between [a/Fe] = 0.4 and [o/Fe] = stars with l^j j = 4000 K, 
[Fe/H] = solar, log(g) = 3.0 from the models of Coelho et al. (2005). 



tation of [a/FcJ-cnhanced SPs (whose predictions do not 
extend far enough to test this). 

5. DISCUSSION 

The above analysis has important implications for the 
modeling and interpretations of integrated stellar pop- 
ulations. The full population synthesis fitting of NGC 
628 and UGC 2124 indicates the presence of a young 
component (~ 1 Gyr, where MarQ5 and BC03 model pre- 
dictions differ significantly) contributing ~ 50% to the 
y-band flux. In SPS models which account for a detailed 
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Fig. 6. — Same as in Figurc|4l but for NIR colors. 



treatment of all stages of the TP-AGB stellar evolution 
phase, this young population would contribute a signif- 
icant amount of NIR flux, resulting in very red optical- 
NIR colors, a mild reddening in J — H, and a bluing of 
the H — K and J—K colors. All of these predicted trends 
were confirmed with the observed colors, thus providing 
further support for the presence of the young component, 
in addition to strongly favoring the contribution of the 
TP-AGB phase used in the Mar05 SPS models (as weU 
as the Marigo & Girardi (2007) models which agree well 
with the Mar05 predictions). 

Such anomalous colors have been reported before. For 
a sample of 5800 galaxies from the SDSS, Eminian et al. 
(2008) compared optical and NIR colors with quantities 
derived from spectra (star formation rate (SFR), age, Z , 
and dust attenuation). They found that galaxies with 
higher SFRs, while having bluer optical colors (g—r), also 
had redder NIR colors {H — K), which they interpreted 
as being in qualitative agreement with the dominance in 
NIR light of the TP-AGB phase after a burst of SF. 



At higher redshift, Elhs, Abraham, & Dickinson (2001) 
observed the colors of bulges out to z 1 and found 
that, while many of the bulges had blue optical colors, 
consistent with recent SF, at z > 0.5, the same bulges had 
red observed J — H colors. At these higher redshifts, this 
translates roughly into rest-frame i — J. Indeed, while the 
Mar05 models only show a mild reddening compared to 
those of BG03 of the J—H color (by '^0.1 mag for a solar 
1 Gyr SSP), the i — J color is much redder (by ~ 0.7mag), 
thus explaining why the effect is only observed at higher- 
z. While Ellis et al. attributed the red NIR colors to a 
predominantly "short burst" mode of bulge building at 
high-z, the red NIR colors could also be explained by the 
presence of young SPs in the active TP-AGB phase. 

These results indicate that caution must be taken when 
interpreting the optical-NIR colors of integrated stellar 
populations. Red colors have typically been interpreted 
as being due to some combination of old age, high Z, 
and significant dust extinction. However, we have shown 
here that some of the redness could actually be due to 
a significant contribution to the SED from a young SP 
with significant NIR fiux excess from the TP-AGB evo- 
lutionary phase. Further confirmation of this result and 
its implications awaits larger samples and analysis with 
forthcoming higher-resolution implementations of the SP 
models which include proper treatment of the TP-AGB 
stellar evolutionary phase. 
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